Schrodinger Cat States of a Nanomechanical Resonator 
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We present a scheme of generating large-amplitude Schrodinger cat states and entanglement in 
a coupled system of nanomechanical resonator and single Cooper pair box (SCPB), without being 
limited by the magnitude of the coupling. It is shown that the entanglement between the resonator 
and the SCPB can be detected by a spectroscopic method. 
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The fabrication and probing of ultra-small nanonie- 
chanical resonators with secular frequencies of GHz and 
quality factors approaching 10^ have been achieved in re- 
cent experiments These resonators are promising sys- 
tems for demonstrating the quantum mechanical nature 
of the mechanical degrees of freedom 2] • Potential appli- 
cations include detection of weak forces, precision mea- 
surement, and quantum information processing 'J 0, |^. 
One crucial step in studying the nanomechanical res- 
onators will be the quantum engineering and the de- 
tection of the mechanical modes. This can be achieved 
by connecting the resonators with solid-state electronic 
devices 0, 0, 0; fo'^ example coupHng a resonator with 
a single electron transistor (SET) via electrostatic in- 
teraction. The SET measures the fiexural oscillation of 
the resonator with an accuracy approaching the quan- 
tum limit 01 . Cooling of a resonator to its ground state 
has been proposed by quantum feedback control via a 
SET0 and by side band cooling via a quantum dot 9] or 
aSCPB[i3. 

The resonator modes can be treated as underdamped 
harmonic oscillators with the damping described by the 
finite quality factor Q. Connecting a resonator with a 
quantum two level system forms a spin-oscillator model 
which has been intensively studied in quantum optics, 
especially in ion trap quantum computing I IjJ. Hence 
the techniques of manipulating the motional state of a 
trapped ion by laser control of its internal mode can be 
applied to studying the nanomechanical resonators lllj. 
In Ref . Q , the capacitive coupling between a nanome- 
chanical resonator and a SCPB was studied, where the 
SCPB can be treated as a two level system - the su- 
perconducting charge qubit - by adjusting the param- 
eters and the gate voltage It was shown that en- 
tanglement between the resonator and the qubit can be 
generated and be detected by interferometry when the 
coupling is stronger than the energy of the resonator. 
In this paper, we show that large-amplitude Schrodinger 
cat states and entanglement can be generated in the cou- 
pled resonator and SCPB system by parametric pump- 
ing of the SCPB when the magnitude of coupling is weak 
due to the geometry of the charge island and the dis- 
tance between the charge island and the resonator 0. 
Given the large amplitude of the generated cat states. 



the entanglement between the resonator and qubit can 
be observed by a spectroscopic measurement which se- 
lectively flips the charge qubit depending on the state of 
the resonator. When the scheme is generalized to two or 
more nanomechanical resonators, it generates maximal 
entanglement between these resonators, which is a key 
element in continuous variable quantum computing[l3j. 
In precision measurement, the cat state of N resonators 
can increase the sensitivity to weak forces by a factor of 
VTV 4]. The effect of environmental fluctuations includ- 
ing the mechanical noise and the charge noise around the 
charge island is analyzed. Furthermore, as this scheme 
involves the generic system of one spin and one oscilla- 
tor, it can be tailored for other applications such as single 
spin detection by a cantilever jl^. 
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FIG. 1: Left: a resonator couples with a SCPB. Right: the 
trajectory of the resonator starts from the origin x — (thin 
dotted line) and the state | j) of the charge qubit; the qubit 
is flipped every half period of the resonator n/uio- 

The coupled system of a nanomechanical resonator 
and a SCPB is shown in Fig. ^ with the resonator un- 
dergoing fiexural vibration. The fiexural mode is de- 
scribed by the Hamiltonian iJ.,„ = huiQa^a where ujq is 
the frequency of the mode and a} (a) is the raising (low- 
ering) operator of the mode. The resonator is biased 
at a voltage Vx{t) and couples to the SCPB through a 
capacitance Cx{x) = C^oi^ + x/cLq) where C^o is the 
static capacitance, dg the static distance between the 
two, and x — (5xo(a -\- a)) the displacement of the fiex- 
ural mode with 5xq — ^fijlmLOQ the quantum width 
of the resonator. The SCPB is a superconducting is- 
land connected with Joscphson junctions and is con- 
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trolled by the gate voltage Vg{t) through the gate ca- 
pacitance Cg. When CgVg + CxVx = (2n + l)e + 2e5n 
with n integer and 5n ^ 1, the SCPB can be treated 
as an effective quantum two level system - the super- 
conducting charge qubit^^l ~ described by the Hamil- 
tonian Hq = 4EJnaz + ^^o^ with = /2C^ 
the charging energy, Ce the total capacitance connected 
with the charge island, and Ej{t) the Josephson energy 
with a small modulation around the static Josephson 
energy Ejq. Here ax,z are Pauli operators for the two 
level system. To the lowest order, the coupling between 
the resonator and the SCPB is —A {t) (a -I- a)) Gz with 
A(t) = 4i?c(14(t)Ca;o/2e){feo/do), and its magnitude can 
be limited by Cx^ and the small ratio of feo/do- In our 
scheme, the bias of the resonator is 14 (i) = cos {ujact)-, 
an ac voltage with amplitude T4o and frequency LOac\ the 
gate voltage is Vg {t) = Vdc + VgQ cos {ujact) , including a 
dc voltage Vdc and an ac voltage with amplitude Vgo and 
frequency uJac- We let CgVdc = (2n -I- l)e so that the 
charge qubit operates at the degenerate point (15|. In ex- 
periments, it was shown that the decoherence time of the 
charge qubit at the degenerate point can be as long as 
microseconds 0. Here uiac = Ejq/U is in resonance with 
the energy of the charge qubit at the degenerate point. In 
the following we study the system in the rotating frame of 
^^Ux- In this frame the Hamiltonian after the rotating 
wave approximation can be derived as 



Hrot — Tlloqo) a 



Ao 



(1) 



where Aq = AEc{VxQCxo/'2e){Sxo/do) is the coupling in 
the rotating frame, = 8Ec{CgVgo + CxoVxo)/'2e < Ec, 
and e_L — Ej{t) — Ejq is the small modulation of the 
Josephson energy. Here, the dynamics of the resonator 
is that of a shifted harmonic oscillator with the Hamilto- 
nian DH,nD^ with the displacement operator 



D = cxp — (a — a 



2hLU(] 



Typical parameters [lOj are Ejo ~ 10 GHz, Ec ~ 50 GHz, 
Cxo ~ 20 aF, and loq « 100 MHz. With Vxq ~ 1 V, the 
coupling is Aq « 20 MHz. 

Below we show that by pumping the charge qubit with 
stroboscopic pulses, Schrodinger cat states with large am- 
plitude can be generated in this system. In an ideal sit- 
uation, we consider 5-function pulse sequence 



e± (t) ~ ttS (riTo) , withyi > 1, integer 



(2) 



where each pulse is a transformation —iax that flips the 
charge qubit after every half period of the resonator mode 
To = tt/wq. Here = 0. The (5- function approximation 
is vahd when e± 3> hujo,^o- Let Ui = e"*^"^"'"^" |e^=o 
be the free evolution of the system between the pulses: 
Ui = De~^'^'^^ ""DK At times utq after the nth pulse, 
the unitary transformation on the system is U (nro) = 



{—idxUiY'. With the relations: UxDcTx = -D^ and 
giTrata^jg-JTrata = ^j'f, we derive 

{D^f"^, n G even 

UinTo)^{ . (3) 
cr^e-"" n G odd 

where the overall phase factors are omitted. This trans- 
formation generates in the wave function of the resonator 
a displacement of Ax = —Sxq {2n\Q<j ^ I Tiloq) when n is 
even, and an opposite displace when n is odd. With the 
initial state = (co| T)+ci| i)) / dx\x)ip{x), where the 
states I ti i) are eigenstates of the charge qubit in the az 
basis and (f{x) is the wave function of the resonator, after 
even number of flips n, the state is 

Colt) / dx\x) ip{x + Ax) + ci\i) / dx\x)ip{x-Ax) (4) 



where the state of the resonator is shifted according to 
the state of the charge qubit. 

Assume an initial state of |^o) = -^(1 T) + I i))|0) 
where |0) is the ground state of the resonator and Ao = 
at i < 0. Following Eq. after even number of pulses 
n, the state ^(| t )| - 2naa) + | [ )\2naQ)) is gener- 
ated where ao = Xo/2huJo is the dimensionless coupling. 
The state \a) denotes a coherent state of the resonator 
mode with an amplitude a. When 2?iao 3> 1, maximal 
entanglement is generated between the resonator and the 
charge qubit. An intuitive way to describe the process is 
to consider a classical particle with two spin components 
in an harmonic potential, where the potential is shifted 
from the origin to the left (right) at spin down (up) when 
the coupling is on. The spin is subject to kicks every half 
period of the oscillator, as shown in Fig. At i < 0, 
the oscillator is at the origin in its ground state with no 
coupling. At i > 0, the coupling is turned on and the par- 
ticle starts oscillating with opposite displacements for the 
two spins. Each kick maximally increases the energy of 
the particle and generates coherent states with large am- 
plitude. Writing the generated state in the |±) basis, we 
have i|+)(| - 2nao) + |2nao)) + i|-)(| - 2nao) - \2nao)). 
A measurement on the ax operator of the charge qubit 
with the scheme in ^ij projects the resonator to the state 
-^(1 — 2nao) ± |2nao)) corresponding to the measured 
ax value -I- or — respectively. Note that with the Hamil- 
tonian in Eq. entanglement can be generated at the 
degenerate states | T)|q!o) and | i)| — ao) without the 
pumping processpQ. However, with Ao ^ hujQ, the cou- 
pling only slightly shifts the resonator state: ao ^ 1 and 
the resonator is only weakly entangled with the charge 
qubit. With the pumping process, a shift of the wave 
function much larger than ao can be achieved. 

This scheme can be generalized to multiple resonators 
and (or) charge qubits to generate entanglement between 
the resonator modes. When two resonators couple with 
one charge qubit with the coupling J^^i^i + 
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the state ^(| - q;i)i| - a2)2 ± |Q;i)i|a2)2) can be gener- 
ated, where the index i = 1, 2 labels the two resonators 
and ai = nXi/huJo. Such states are maximally entangled 
states of the resonator modes and are crucial elements 
in continuous variable quantum computingfisj. 

The entanglement and coherence between the res- 
onator and the charge qubit can be detected by a spec- 
troscopic method for resonator states of large amplitude 
with 2nao ^ 1- During the detection, choose a static 
bias Cz > of the charge qubit and modulate the Joseph- 
son energy by cos (oJdt) Ox with a frequency LOd for a 
duration of i:/uJd- Here instead of the i5-function pulses 
in Eq. 12 Cd has the same order of magnitude as SnaoAo 
and e^, while ^ wq and the resonator can be treated 
as static during the detection. This condition is crucial 
in realizing the detection process. The effective Hamilto- 
nian of the SCPB is then: 

r7-|±2nQn) ^ 4naoAo , , , , 

for the resonator states of |±2nao) respectively. Hence 
the energy splitting of the charge qubit depends on the 
state of the resonator with a splitting i5+ = + 4naoAo 
for the resonator state |2nao) and E- — tz — 4naoAo 
for the resonator state |— 2nao). We choose the pulse 
frequency to be fiiOd = (-z — 4naoAo, in resonance with 
the charge qubit in li\ ^""o) _ ^his pulse is then followed 
by a (5-function 7r/2 pulse that transforms | t, i) to | + , — ). 
Applying the pulses to the state -^(| \ )| — 2nao) + 
I I )| 2n ao )), the final state is 

-^h>l - 2nao) + (cj-) -f ctl-f)) \2na^) (6) 

with C| = -isin (7re<j/2ed) {ed/ld) and 

/TTErfN . . (TTed\ SnaoAo 

c I = cos - — — I sm - — 

\2td) \2edJ Cd 



where = y -I- (SnaoAo) . For the resonator state 
|2nao)j the off resonance SnaoAo between cod and 
prevents the charge qubit from flipping. By adjusting 
the bias Sz and the amplitude e^, we can find a regime 
where |cj « 1. Note the states |±) are in the rotat- 
ing frame and in the lab frame the ax eigenstates are 
|±)s = e**^-'*/''|±). A measurement on the ax operator 
of the charge qubit as in |l5 | obtains the probabilities 
of the states |±): p_ = (l-f |cxP)/2 and p+ = \c^\^/2 
respectively. As a first step, the correlation between the 
resonator and charge qubit can be demonstrated by this 
measurement when p_ ~ 1 and p+ ~ 0. When no corre- 
lation exists between the resonator and the charge qubit, 
~ and p± = 1/2. 

By measuring the charge qubit, it can also be shown 
that the states | T)l — 2nao) and | J,)|2nao) are in coher- 
ent superposition. This measurement starts by applying 



a (5-function ■k/2 pulse to the state -^{\ t ) I — 2n,Q!o ) -f 
11)1 2nao ) ) , followed by the pulses in Eq. |5J for n time. 
The state becomes + \-)\tp_)) with !?/;+) = 

I - 4nao) + 2|0) - |4nao) and \^-) = \ - Anao) + |4nao). 
The probabilities of the states |±) are hence p+ = 3/4 
and p- — 1/4. Without the coherence, i.e. that of a 
mixed state of | t)l ~ 2nao) and | J,)|2nao), the probabil- 
ities are p± = 1/2 respectively. Hence measurement of 
the ax operator probes the coherence of the system. 




FIG. 2: The main plot: the fidelity of the amplification versus 
e± for n = 4, 8, 12 pulses from top to bottom. Inset: the 
probability p_ versus at — 4.0ci;o (solid line) and tz = 
3.2LU0 (dotted line). Here Ao = 20 GHz. 

Ideal situations are assumed in the above discussions 
with well separated energy scales; ej_ 3> ujq during the 
amplification and Ez, e^, SnaoAo ^ ujq during the de- 
tection. In practice, the frequency of the resonator is 
around 100 MHz; while the amplitudes of the pulses e±, 
ed are upper bounded by the Josephson energy Ejq of 
the charge qubit which is typically below 20 GHz. The 
dynamics of the resonator may have important effects 
on our scheme. Below, we numerically simulate the dy- 
namics of the amplification and detection with the above 
parameters and in the coordinate space of the resonator. 

Let the wave function at time t be \ip{t)) — 
J2s I dx\s)\x)(ps{x,t) where s is the state of the charge 
qubit and (/?s(a;,t) is the wave function of the res- 
onator, with the initial state discussed above. We cal- 
culate the fidelity of the amplification process: f{e±) = 
lii'idit) |V'(*))I : where |V'id) is the ideal wave function 
generated by the pulses in Eq. (0). In FigEl the fidelity 
is plotted versus the amplitude of the pulse £± . It can be 
seen that at ej^ = IOujq, the fidelity can be very low with 
f = 0.7 after n = 12 pulses; but the fidelity increases 
rapidly with increasing ej^. At ej^ = 60a;o, correspond- 
ing to e_L = 6 GHz, the fidehty is / > 0.99 after n = 12 
pulses. In the detection process, the dynamics of the 
resonator affects the probability We simulate the 
detection process at various static bias Cz and with in 
the range of 0.5 wo — 10.5 cjo. Here SnaoAo ~ 1.9 wo after 
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n = 12 pulses. In the inset of Fig|21 P- is plotted versus 
Ed at £2 = 3.2 cjQ and ez = 4.0(^0. When ^ SnaoAo, 
the charge qubit flips for both the states | ± 2nao) and 
p_ ^ 0.5. When ^ SnaoAo, the off resonance strongly 
affects For = 4.0 ljo, a maxumum of p- appears 
at Erf ~ 1.9ciJo with p_ = 0.80, very different from the 
probability without the correlation. For = 3.2 o^o, p^ 
oscillates with Cd due to the evolution of the resonator. 
This shows that the correlation between the resonator 
and the charge qubit can be detected even at high res- 
onator frequency. 

With a couphng of 20MHz[l3|, n > 10 flips are re- 
quired to have 2nao > 1, a duration around 50nsec. It is 
crucial to have a decoherence time longer than this du- 
ration to successfully generate the entanglement. During 
the amplification process, the charge qubit operates at 
the degenerate point where the charge noise, dominated 
by the low frequency charge fluctuations, causes a deco- 
herence time of microseconds [l5j| . In the rotating frame 
of Eq. 1^, this can be explained by a spectral shift: the 
noise spectrum in this frame is S'^ (oj zt Ejo/Ti) with a 
shift of Ejo from the noise spectrum in the lab frame 
S'^ (uj). The low frequency noise is hence screened by 
the Josephson energy. The mechanical noise of the res- 
onator can be described by the quality factor Q. At 
temperature T = 20mK with Q = W^, the dissipa- 
tion rate is ksT/Q = 50 KHz. The decoherence rate 
is rj^^ « {2naQ)^kBT/Q, which at 2nao ~ 5 gives 
'''dec ^ ^ MHz and limits the amplification process. Mean- 
while, in a situation where the phase coherence between 
the states | T)l — 2nao) and | l)\2nao) is not required 
as in the spin detection0|, the amplitude of the gen- 
erated coherent states can be bounded by the quality 
factor. Assume after flippings the amplification sat- 
urates. This means that starting from the coherent 
state I — 2nsaQ) with the charge qubit at | t)? after a 
time of tt/wq, the coherent state is |2nsQ!o). The ini- 
tial elastic energy is Tiwo^o + 1)^7 while the final 
elastic energy is hoJoaQ (2n<, — 1)^ with an energy loss of 
SE = SushujoaQ. The loss is caused by the dissipation: 
SE — AirnlhujQa^/ Q, from which the saturation limit: 
lis — 2QI'K can be derived. 

This scheme involves a generic model of one harmonic 
oscillator and one quantum two level system (spin) cou- 
pling via linear interaction xoz, and hence can be gen- 
eralized to other spin-oscillator systems. One example 
is the single spin detection by magnetic resonance force 
microscopy (MRFM)[l3] where the spins near a surface 
interact with the magnetic particle attached to a can- 
tilever and affect the vibration of the cantilever. By ob- 
serving the frequency or amplitude of the vibration with 
optical interferometry, the distribution of the spins can 
be detected. In experiments, the resolution of MRFM 
has been improved towards the single spin levelHi. In 
the conventional scheme, the cyclic adiabatic inversion 



(CAI) method []J is applied where the spins are driven 
by continuous microwave with periodic modulation of the 
phase of the microwave. Our scheme provides an alter- 
native to this approach. By applying parametric pulses 
to flip the spins every half period of the cantilever, co- 
herent states of the cantilever with large amplitude can 
be generated within a short time even at weak coupling. 
Using the same notations as that in Eq. and assum- 
ing a total local spin of m/2 near the tip, after n spin 
flips in Eq. the coherent states are | ± 2nmao). This 
shows a resolution of 6a/ 6m = 2nao. When 2nao S> 1 
and Us > n, single spin resolution can be achieved. This 
requires Q > tt/Auq of the cantilever. 

We studied a scheme of generating and detecting 
Schrodinger cat state in the coupled resonator and SCPB 
system. Compared with previous works j^i2|, large am- 
plitude coherent states can be obtained at much smaller 
coupling than the energy of the resonator. The scheme 
provides a practical way of investigating the quantum 
properties of the nanomechanical resonators. 
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